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From bcc to fee: interplay between oscillating long-range and repulsive short-range 

forces 
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This paper supplements and partly extends an earlier publication, Phys. Rev. Lett. 95, 265501 
(2005). In d -dimensional continuous space we describe the infinite volume ground state configu- 
rations (GSCs) of pair interactions (fi and (p + ^p, where ip is the inverse Fourier transform of a 
nonnegative function vanishing outside the sphere of radius Ko, and t/j is any nonnegative finite- 
range interaction of range ro < jd/Ko, where 73 = y/6-jT. In three dimensions the decay of ip can 
be as slow as ~ r-~^, and an interaction of asymptotic form ~ cos(A'or + 7r/2)/r'^ is among the 
examples. At a dimension-dependent density pd the ground state of <^ is a unique Bravais lattice, 
and for higher densities it is continuously degenerate: any union of Bravais lattices whose reciprocal 
lattice vectors are not shorter than Kq is a GSC. Adding decreases the ground state degeneracy 
which, nonetheless, remains continuous in the open interval {pd,p'd), where p'd is the close-packing 
density of hard balls of diameter tq. The ground state is unique at both ends of the interval. In 
three dimensions this unique GSC is the bcc lattice at pa and the fee lattice at pjj — v^/j'o. 

PACS: 61. 50. Ah, 02.30.Nw, 61.50.Lt, 64.70.Dv 



I. INTRODUCTION 

In an earlier Letter [1] we described the infinite volume 
ground state configurations (GSCs) of a class of classi- 
cal particle interactions in d > 1 dimensional continu- 
ous space. These pair interactions have a nonnegative 
Fourier transform vanishing above some finite wave num- 
ber, Kq. We proved that at a threshold density pd cx 
there is a unique periodic GSC (the basic-centered cu- 
bic lattice in three dimensions), and above pd the set 
of GSCs is continuously degenerate and contains peri- 
odic and aperiodic configurations. While this was prob- 
ably the first result providing specific examples in three 
(and higher) dimensions, important rigorous work pre- 
ceded it in lower dimensions; see, for instance, Kunz [2] 
on the one-dimensional one-component plasma, Ventevo- 
gel, Nijboer and Ruijgrok [3] and Radin [4] on ground 
states in one dimension and Theil's recent proof of ground 
state crystallization in two dimensions [5]. Although we 
will be concerned only with ground state ordering, let 
us note that rigorous results on phase transitions or or- 
dering in a continuum at positive temperatures do not 
abound, and all are about more or less contrived model 
systems. The one-component plasma in one dimension is 
ordered at all temperatures [2] ; Ruelle [6] proved segrega- 
tion in a two-component system (the Widom-Rowlinson 
model); Lebowitz, Mazel and Presutti [7] proved vapor- 
liquid transition for particles with two-body attractive 
and four-body repulsive interactions near the mean-field 
limit; in a one-dimensional model with an unstable inter- 
action low-temperature freezing into an ordered configu- 
ration was shown by the present author [8] ; and recently 



Bowen et al. proved fluid-solid phase transition in a two- 
dimensional system of decorated hard hexagons [9]. 

In this paper we supplement and partly extend the re- 
sults of [1]. Apart from recalling the definitions, we do 
not repeat what is written there. The main part of the 
theorem of [1] will be stated and proven in a new, simpler 
form, emphasizing the nice algebraic structure of the set 
of GSCs. Moreover, the results will be extended to inter- 
actions of a non-integrable decay. The proof in [1] was 
based on the Poisson summation formula. This formula 
is widely used in physics; one of its earliest and most fa- 
mous applications was the calculation of the Madelung 
constant of ionic crystals by Ewald [10]. The formula in- 
volves at least one infinite summation, and neither the 
convergence of the infinite sum(s) nor the equality of the 
two sides is guaranteed. Although the results of [1] were 
already formally valid to the larger class of interactions, 
we stated them only for a restricted class, those of the 
strongly tempered interactions (see later), because no ar- 
gument supporting the applicability of the Poisson for- 
mula to functions of a non-integrable decay was given 
in that paper. The extension of this formula constitutes 
an active field of research in mathematics, see e.g. [12], 
but is not our main concern here. Therefore, without 
looking for the most general formulation, we propose an 
extension just suitable for our purposes. As a matter of 
fact, the extension involves also the notion of a ground 
state configuration. The definition of a GSC is based on 
infinite sums that are absolutely convergent for strongly 
tempered interactions, but only conditionally convergent 
for interactions of a non-integrable decay, and the way 
they converge has to be specified. 



2 



Another, gratuitous, extension, mentioned but not ex- 
ploited in [1], will be obtained by modifying the short- 
range behavior of the interaction. The inverse Fourier 
transform of an integrable function is bounded and 
continuous — this is our case. Such bounded functions 
play a role as soft effective interactions in polymer physics 
[13], but not in traditional solid state physics where Pauli 
exclusion gives rise to a practically infinite repulsion at 
overlaps of atoms. Imagine, however, that an infinite con- 
figuration X was shown to be a GSC of 9?. Then X will be 
a GSC of all interactions ^p + tp, where il> is non- negative 
and vanishes at and above the nearest-neighbor distance 
oi X: tp does not contribute to the specific energy of X, 
and can only increase the energy of any perturbation of 
X. Reversing the argument, we may start with (f + tp, 
where tp is of bounded support, non-negative and may 
contain a hard core or diverge at the origin as fast as 
we wish. Then (f + tp has common GSCs with (f if the 
support of tp is small enough and the density is not too 
high. 

The physical importance of the above two extensions is 
that with them we obtain the GSCs of interactions whose 
asymptotic form is ^ cos Kor/r^, as that of the RKKY 
interaction, but which can be arbitrarily strongly repul- 
sive at small distances. As noted also by Likos [14] , such 
interactions can model those between ions in metals and 
be relevant in the explanation of the crystal structure of 
certain metals. However, further study is necessary be- 
fore any conclusion could be drawn about this question. 
As an immediate gain, we will find that at some den- 
sity pjj > p3 the unique GSC oi (p + ^p is the fee lattice, 
while at ps it is the bcc lattice. This transition from bcc 
to fee with an increasing density is the consequence of 
an interplay between a long-range oscillating interaction 
(which is short-range in Fourier space) and a short-range 
positive pair potential. 

The following section is the central part of the paper. 
After introducing the necessary definitions we enounce a 
theorem in a rather compact form, and then expand its 
content in a series of remarks. The Poisson summation 
formula is presented here as a lemma. In Section III we 
prove the lemma, an auxiliary statement about Bravais 
lattices, and the theorem. This section also contains the 
proof of a general assertion about the non-existence of 
metastable ground states. The paper is closed with a 
brief Summary. 



II. DEFINITIONS, NOTATIONS, RESULTS 

We consider a system of identical classical particles in 
M'', that interact through translation invariant symmet- 
ric pair interactions, (p{r — r') = ip{r' — r). Rotation 
invariancc is not supposed. An A^-particle configuration 
(iV < 00) is a sequence (ri, . . . , rjv) of points of M'' 
and will be denoted by B (referring always to a Bravais 
lattice), R, X and Y. While the order of the points is 
unimportant, two or more particles may coincide in a 



point, resulting = • • • = r^^^^. Such a coincidence can 
occur if (p{0) is finite, and it indeed occurs in certain 
GSCs to be described below. Throughout the paper, the 
notation (p will be reserved to bounded interactions; un- 
bounded interactions, such as those diverging at the ori- 
gin or including a hard core, will be composed as p + tp. 
The number of points in R will be denoted by Nr. The 
energy of a finite configuration R is 

u(R) = l E ^i^-^')- (1) 

r,r'G-R,r^r' 

Let i? be a finite and X be an arbitrary configuration. 
The interaction energy of R and X is 

/(i?,X) = ^/(r,X) = ^^V'(r-x), (2) 

and the energy of R in the field of X is 

U{R\X) = U{R)+I{R,X). (3) 

If X is an infinite configuration, the infinite sum in (2) 
has to be convergent. This imposes conditions on both X 
and 93, and the stronger the condition on X, the weaker 
it can be on p. For example, one may ask I{R, X) to be 
finite for every X that is locally uniformly finite, meaning 
the existence of an integer mx such that the number of 
particles in a unit cube everywhere in stays below mx ■ 
This was our choice in [1]: the corresponding condition 
on the interaction is strong temperedness which for a 
bounded ip reads 

E i^wi < °° (4) 

for any locally uniformly finite X. 

Definition. — Given a real /i, X is a ground state config- 
uration of p> for chemical potential /x (a /uGSC) if for any 
bounded domain A and any configuration R 

uiRn^\x\^)-^iNB.nA > u{xnA\x\A)~^iNxnA (5) 

where X OA and X\A are parts of X inside and outside 
A, respectively. X is a ground state configuration (GSC) 
if (5) holds true for every R such that N^da = NxnA- 

A seemingly more general, but actually equivalent def- 
inition is as follows. X is a /iGSC (respectively, X is 
a GSC) if for any finite part Xf of X and any finite R 
(respectively, any R such that = Nxf ) 

U{R\X \Xf)-fiNR> UiXf\X \Xf)- fiNx, • (6) 

If p} is strongly tempered, we can — at least in princi- 
ple — test any locally uniformly finite X to be, or not, 

a GSC according to (6). Ground states of interactions 
that violate condition (4), as those between ions in met- 
als mediated by the Priedel oscillation of the conduction 
electrons, can be defined only within a more restricted 
set of configurations. Intuitively, ground states cannot 
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be arbitrary sets of points, they are arrangements with 
some good averaging (ergodic) property. Specifically, we 
shall look for them only among periodic configurations 
and their unions. Simultaneously, the infinite sums ap- 
pearing in (3) will be suitably interpreted. 

A Bravais (direct) lattice B = {X^^^i fiaaajn e Z**} 
is regarded as an infinite configuration. Here 

are 

linearly independent vectors and n = {ni,...,nj_) is 
a d-dimensional integer. The dual (reciprocal) of B 
is the Bravais lattice B* = {^naba|n e Z"^} where 
• b/3 = 27r(5a/3. The nearest neighbor distances in B 
and B* arc denoted by rs and qs', respectively. The lat- 
ter is related to the density of B via p{B) = C"b"{<1b')''', 
where C"b" is determined by the aspect ratios and an- 
gles of the primitive cell of B. [Notational remark: B 
and B* will always refer to specific Bravais lattices as 
given above. "S" refers to the family of all Bravais lat- 
tices of the type of B, characterized by dimensionless 
quantities. "_B" may take on the 'value' bee, fee, simple 
cubic, and so on.] We shall look for GSCs of the form 
X = Lij^i{Bj + Yj) where Bj are Bravais lattices and 
Bj + Yj is Bj shifted by the vector yj . We shall refer to 
configurations of this form as unions of periodic configu- 
rations. If J = 1, X is a Bravais lattice. If Bj = B for 
each j then X is periodic. If at least two diff'erent Bravais 
lattices are involved in the union then X is either peri- 
odic or aperiodic, and may contain overlapping points 
that are to be counted with repetition. The density of X 
is p(X) = ELp(S,). Now 



7(r,X) = ^7(r-y,-,i?,), 

i=i 



(7) 



so the sum to be interpreted is 



I{r,B)= ^ ^{r + -R) = 

RGB 



(p + ^^nc,aLaj . (8) 



The interaction ip will be defined as the inverse Fourier 
transform of a function ip ^ (M.'^') that vanishes outside 
the ball of radius Kq: (/?(r) = {2tt)-'^ J^^j^^ ^(k)e»''"" dk. 
Then (p is continuous and all its derivatives exist and are 
also continuous functions decaying at infinity; in fact, 
(p{r) is an entire function of r [11]. By definition, 

^ ^(r + R) = hm ^ e-^l-+^lV(r + R) (9) 
Res Res 

provided that the limit exists. In our case a weaker, 
e.g. exponential, tempering would suffice and give the 
same result; the Gaussian tempering is more convenient 
to work with in an arbitrary dimension. Note that the 

sum on the right-hand side is absolutely convergent for 
any e > 0. It is easily seen that whenever the sum in (8) 
is absolutely convergent, Eq. (9) yields the same result. 
This is the case of all the examples given in [1]. In the 
absence of absolute convergence, the sum (8) can still be 



conditionally convergent; e.g. with some mild additional 
assumption on p one can show that 



/(r, B) = lim 

JVi,...,Ar<i^oo 



E 



neZ<*,|na|<iV„ 



(10) 

exists and agrees with the result suggested by the Poisson 
summation formula. However, the proof of this formula 
is simpler with the definition (9). 

THEOREM. Let (p € ^^(M'*) be a real function with 
the following properties: 

(1) p is continuous at the origin, 

(2) p{~-k) = p{-k), 

(3) > and 

(4) there is some Kq such that ^^(k) = for |k| > Kq. 
(i) Define p{r) = {2tt)-'^ J <^(k)e*'^ "' dk. Choose Bra- 
vais lattices Bi,. . . ,Bj such that each qb* > Kq, where 
equality is allowed only if p is continuous at |k| = Kq. 
Then X = u/^i(Bj + y^) is a GSC of p for arbi- 
trary translations yj and it is also a /xGSC for /i = 
p{X)ip{0) — ^p{0). The energy per unit volume of X 
is e{X) = e{p{X)) where 



4p) = Ip[p0{o) - V(0)] 



(11) 



is the minimum of the energy density among unions of 
periodic configurations of density p. GSCs of the above 
properties exist in a semi-infinite density interval [pd, oo). 

(ii) If p is strongly tempered and X is locally uniformly 
finite with existing p{X) > pd and e{X) > e{p{X)), then 
X is not a GSC. If (p is not strongly tempered but the 
limit (10) on Bravais lattices exists, then any union X 
of periodic configurations with p{X) > pd and e{X) > 
e{p{X)) is not a GSC. 

(iii) Let ro < ^d/Ko, where 71 = 27r, 72 = 47r/\/3 
and 73 — VGtt, and let t/j be a real function such that 
ip{r) G [0,00] and ^(r) = for r > ro. If vb > ro and 
te* > Kq, then B is a GSC and a /uGSC of (p + tp with p, 
and e{B) given above, not depending on ip. GSCs of the 
above properties exist in a density interval [pd, p^]- 

Remarks. — 1. Compared with the theorem of [1], the 
conditions on ip are formulated uniquely via p, and are 
considerably weaker. For instance, p or its derivative can 
be discontinuous at Kq. Here are two examples in three 
dimensions: (,S(k) = 1 for A: < iCo yields (fc = |k|, r = |r|) 

p{r) = -{Ko/2Tr^)cosKor/r'^ + {l/2TT'^)smKor/r^ , 



while with (^(k) = 1 — k/Ko fov k < Kq we obtain 

, , cos(/v()/- + 7r/2) 1 — cos/vi)/- 

'^(0 = ^Z23 + 



TT^i^or'' 



(12) 



(13) 



Both arc conditionally summable on Bravais lattices [if 
B* has no point on the sphere |K| = Kq, in the case of 
(12)], as defined in Eq. (10). 

2. We proved in [1] that the condition qs* > Kq can be 
satisfied only if p> pd, a dimension-dependent threshold 
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density at which qb* = Kq for a unique Bravais lattice 
B, and this is the unique periodic GSC; in particular, 
P3 = Xq/Sv^tt'^ and the lattice is the bcc one. The 
above form of the theorem shows that for Pd < P < '^Pd 
no union is available, only Bravais lattices can be GSCs. 
In [1] we gave also the densities of some Bravais lattices 
B at which qb* = Kq. Recalling these values, 

4\/2 

Pbcc = P3 < Pfcc = P3 = 1.089p3 

< Psh = P3 < Psc = \/2p3 (14) 

(sh = simple hexagonal with c/a = \/3/2, sc = sim- 
ple cubic), one can sec that all the high-symmetry Bra- 
vais lattices appear as GSCs between ps and 2p3. Also, 
if p{Z) denotes the density of a metal of valency Z 
then, in the free-electron approximation and supposing 
a spherical Fermi surface of radius kp — Kq/2, p{Z) = 
{V2n/3Z)p3 = (1.481/Z)p3 which for Z = 1 is in this 
interval. In general, in the interval npd < p < {n + l)pd 
the ground state configurations are unions of at most n 
Bravais lattices, each of density > pd- Thus, the simplest 
aperiodic GSCs, unions of two incommensurate Bravais 
lattices, appear only if p > 2p4. For example, in 3 di- 
mensions at 2/93 they arc the unions of two bcc lattices 
rotated and possibly shifted with respect to each other. 

3. The family of all the GSCs of above the den- 
sity Pd is closed on unions. This is obvious from the 
present formulation, because the union of two GSCs of 
the form given in the theorem is a configuration of the 
same form, so it is necessarily also a GSC. Recall from 
[1] that a periodic configuration X is called B-periodic if 
X = Uj^i{B + yj) and B is chosen so as to minimize J. 
If Xm are B^-periodic configurations then LlX^n is peri- 
odic if and only if B = (iBm is a rf-dimensional Bravais 
lattice. Because B C Bm, B* D and qb* < Qb-^- It 
follows that a B-periodic configuration X = [Jj^i{B+yj) 
can be a GSC even if g^. < Kq, provided that it can be 
written also as X = U^^i(Sj +yj), where gs* > Ko for 
j = 1, . . . , J'. This means that on average the Bjs are 
denser than B and thus J' < J. Logically, if we permit 
different Bravais lattices to occur in the union forming 
a periodic configuration, the number of components may 
be decreased. 

4. A pGSC is, by definition, also a GSC because it sat- 
isfies a stronger condition. Therefore, in the theorem it 
would have been enough to say that X or B is a p,GSC. 
We wanted to emphasize that the theorem strengthens 
that of [1] by stating that the opposite is also true: a 
GSC is always a pGSG, even if (p{0) = and thus ip is 
not supcrstable. In the latter case p, = —^(^(0), indepen- 
dently of the density of the ground state configuration. 

5. Assertion (ii) of the Theorem extends to a larger 
class of configurations and pair potentials earlier results 
by Sewell [16] and Sinai [17] on the absence of metasta- 
bility for strongly tempered interactions. Following the 
usual definition [17], we apply the term "ground state 



configuration" as a synonym of a locally stable configu- 
ration. Thus, a GSC could be globally unstable, meaning 
that by some perturbation involving infinitely many par- 
ticles its energy density could be decreased. Such a GSC 
might be called metastable. However, we have to pre- 
cise the kind of infinite perturbations we allow. If the 
particle density is allowed to vary, usually the absolute 
minimum of the energy density is attained at a single 
value of p, and all GSCs of a different density should be 
considered metastable. In this sense, the unique stable 
GSC of an everywhere positive interaction is the vac- 
uum, and for the interactions (p studied in this paper the 
globally stable GSCs are the /xGSCs belonging to ^ = 
(hence, to p = (p{0) /2ip{0), if this value is finite and not 
smaller than pd). We adopt a more restrictive definition 
of metastability, not allowing the density to vary. Then 
the configurations characterized by the theorem as GSCs 
are not metastable because their energy density is the 
attainable minimum for their density, and no other con- 
figuration (within the specified class) can be a metastable 
ground state. 

6. A sufficient condition for (p to be strongly tempered 

is that \(p{r)\ < Cr~'^~^ for r > r' , where C, 77 and 
r' are some positive numbers. In our case this holds, 
for example, if besides conditions (l)-(4), <^ is 3 times 
differentiable, see Eq. (9) of [1]. 

7. Point (iii) of the theorem needs more an explanation 
than a formal proof. When we ask tb > ''0, we limit the 
role of tjj to reducing the degeneracy of the GSCs of (p. 
The largest allowed range of tp, tq — jd/Ko, equals the 
nearest-neighbor distance of the imique GSC of ip at the 
density pd- the uniform chain, the triangular lattice and 
the bcc lattice for d = 1, 2 and 3, respectively. The 
theorem makes no prediction if rg > "/d/Ko, because no 
Bravais lattice satisfies both conditions tb > Jd/Ko and 
Qb' > Kq. This is an obvious consequence of rsqB* = 2n 
in one dimension, and of 

f 47r/V3 (rf = 2) 
7. = maxr5to. = |^^ (d = 3), ^'^^ 

with the maximum attained if B is the triangular lat- 
tice in two, and the bcc (or fee) lattice in three dimen- 
sions. We prove (15) in the next section, and suppose 
henceforth that tqKo < jd- For a given lattice type 
"B", tbQb' = 7"B" = 7"B*", independent of the density. 
This implies that the simultaneous inequalities tb > tq 
and qB* > Kq hold in a (closed) interval of the density, 
whose lower and upper boundaries are implicitly deter- 
mined by gs* = Kq and tb = ro, respectively. If p is in 
this "stability interval" I"b" , then i?isaGSCof(p-l-V', 
with energy density e{B) = e{p). If p is not in L^b" , 
then e{B) > e{p) [provided the strict positivity of for 
k < Kq and of V' for r < ro, that we suppose now], and 
B is surely not a GSC oi p + ip ii p is still in the stability 
interval of some other Bravais lattice B: In this case B 
is a GSC with energy density e{B) = e{p), and due to 
the absence of metastability, B cannot be a GSC. I--b" 
shrinks to a single point if ro = 7"b" /Kq and disappears 
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if ro > 7"B" /Kq. In two dimensions 72 = 7tr, so this may 
happen with any Bravais lattice other than the triangu- 
lar one. Because the triangular lattice is self-dual, at a 
given density it has the largest qs- and the largest rs 
among the two-dimensional Bravais lattices. Therefore 
its stability interval 



[P2,p(r-tr = ro)] = [p2,P2] 



V3K^ 2 



(16) 



contains in its interior the stability intervals of all the 
other Bravais lattices. Figure 1 shows the stability in- 
tervals of the triangular and the square lattices for all 
allowed values of ro. If p falls into In- but outside I"b", 
then B is certainly not a GSC oi (p + ip a.t this density. 
In particular, at the two ends of /tr the only GSC is the 
triangular lattice. The uniqueness at the upper value 
r-Q is new, and is due to ip. 




FIG. 1: Stability of two-dimensional lattices. The stability 
intervals of the triangular and square lattices are obtained as 
horizontal cuts of the respective domains. 

The situation in three dimensions is more complicated, 
because 73 = 7bcc — 7fcc- From the expressions 



rbcc(p) = V3/(4p)i/3, 



.(p) = {V2/pf/^ (17) 



of the nearest-neighbor distances one can see that the 
stability intervals for the bcc and fee lattices are 



^bcc 



8\/27r3 



3^3 
4^ 



^fcc 



6\/37r3 



V2' 



(18) 



the lower boundaries being = pbcc and pfcc, respec- 
tively, cf. Eq. (14). Because at a given density the bcc 
lattice has the largest qb* and the fee the largest rs, the 
stability interval of any Bravais lattice is between ps and 
/O3 = \/2/rQ. According to the argument given above, 
at Pa and P3 the ground state is unique; especially, at 
p — V^/tq the unique GSC is the fee lattice. The in- 
tervals /bcc and /fee only partially overlap, and may not 



overlap at all: if 37r/2^/3 < roKo < VQtt, /bcc and /fee are 
disjoint. For tq < \/6tt/Ko the stability intervals of other 
Bravais lattices fill the gap. At Tq = y/En/KQ, /bcc and 
/fee shrink to a single point, pbcc and pfce, respectively, 
and for densities in between no Bravais lattice satisfies 
both conditions qB' > Kq and r^ > tq. In this interval 
the ground state is probably unique and changes contin- 
uously from bcc to fee as the density increases. In Figure 
2 we present the stability intervals of the bcc and the fee 
lattices for all the allowed values of ro. 




Pfcc P:i' 



FIG. 2: Stability of three-dimensional lattices. The stability 
intervals of the bcc and fee lattices are obtained as horizontal 
cuts of the respective domains. 

The GSCs of (p will be found by applying the following 
extension of the Poisson summation formula. 

LEMMA. Let ip e L^{R'^) be a function of bounded 
support, which is continuous at the origin. Let ip{r) = 
{2Tr)~'^ J (^(k)e''' '' dk. Choose a Bravais lattice B such 
that if is continuous at every K in B* . Then 



lim > ( 

elO ^ 

Res 



-=l'-+i^lV(r + R) = p(S) (^(K)e'K- 



(19) 

implying the existence of the limit. 

Observe that the sum in the right member has only 
a finite number of nonzero terms, hence for any B the 
continuity of is to be checked only in a finite number 
of points. In particular, for any B dense enough the only 
point of B* inside the support of (p is the origin, where 
(p is continuous. This is precisely the fact we shall use in 
the proof of the theorem. 



III. PROOFS 



A. Proof of the lemma 



Let 



5e(k) = (47re)-'^/^e 



(20) 
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First, we show that 



^ e-^l-+«-lV(r + R) = p{B) J2 * ^e)(K)e'^-^ 



RGB 



(21) 

for any e > 0. The argument is essentially the same as 
the one we used to prove the Lemma of [1]. Because 
(p{r) is an entire function of r decaying at infinity, both 
e~^'^ (p{r) and its Fourier transform (p * 6^ are functions 
of rapid decrease [15]. Therefore, the infinite sums on 
both sides of Eq. (21) are absolutely convergent and the 
convergence is uniform in r. So both sums define con- 
tinuous functions, that are periodic with periods € B 
and have the same Fourier coefficients. Indeed, multi- 
plying Eq. (21) by e"*^ "" and integrating by terms over 
the unit cell of volume p{B)^^, on the left-hand side af- 
ter summation we obtain ((p * (5g)(K), which is also the 
trivial result on the right-hand side. Because of the com- 
pleteness of the system {e*^ ''|K G B*} in the Banach 
space of integrable functions on the unit cell of B, the 
two continuous periodic functions coincide everywhere. 

Next, we prove that the integrability of (f and its con- 
tinuity at K imply 



lim((^*5,)(K) = ^(K). 

£—►0 



(22) 



Fix any r] > and write 



((^*^e)(K)= / ^(K-q)5,(q)dq 

Jq<ri 



+ <f>{K- q)5e(q) dq = Je,<„ + Je,>r,- (23) 
Jq>n 

If e is small enough then (5e(q) < 1 for ^ > 77, and in 
Je.yq the integrand can be bounded above by |(^(K — q)|. 
Thus, due to the dominated convergence theorem the 
limit and the integration can be interchanged, resulting 
limg^o Je.yri = 0, because limg^o ^£(q) = for 5 > 77. 
On the other hand, 

Je,<r, = <^(K) / (5£(q)dq 

Jq<ri 

+ I [^(K-q)-^(K)]4(q)dq. (24) 
Now lims^o /g<^ '^^(q) dq = 1 and 

/ [^(K-q)-^(K)]5,(q)dq 



<sup|^(K-q)-^(K)|. 



(25) 



Combining the above equations. 



limJ(<^*5,)(K)-^(K)| <sup|^(K-q)-<^(K)|, (26) 



£^0 



q<ri 



Finally, using the fact that the different sums involving 
i^(K) [but not {ip * 5e)(K)!] are finite, we find 



lim 

£^0 



[(^*5,)(K)-<^(K)]e'^ 

Kes» 

< sup|(^(K-q)-^(K)| 

KGB* 



(27) 



holding for all 77 > 0. Thus, we can conclude that the 
left-hand side is indeed zero. This completes the proof of 
the Lemma. 



B. Proof of Equation (15) 

The proof is based on the well-known fact that I-k jqB* 
is the largest distance between neighboring lattice lines 
[d = 2) or planes [d 3) of B, see e.g. [18]. The 
way we proceed is to consider all the types of Bra- 
vais lattices simultaneously and to select, through a se- 
quence of choices, the vectors aj that define the maxi- 
mizer of rsQE*- Let B be any two-dimensional Bravais 
lattice. Choose ai among the shortest vectors of hence 
ai = rs- The lattice line parallel to ai is a line of largest 
density, therefore the largest distance between neighbor- 
ing lines is measured perpendicular to ai. The other 
primitive vector, a2 is selected among the shortest vec- 
tors not coUinear with ai and making an acute angle a 
with ai. Then a > tt/3 (otherwise ai — a2 would be 
shorter, and should replace a.2), and 27r/gB* = 02 sin a 
(thus, 62 = ?B*)- Therefore 



27rai 
a2 sin a 



(28) 



whose maximum on the condition that a2 > Oi and 
7r/3 < a < 7r/2is 47r/\/3, attained with the choice 
a2 = ai and a = w/S, characteristic to the triangular 
lattice. 

In three dimensions, given B, let P be a lattice plane 
of highest density, containing the origin. Let ai be one 
of the shortest lattice vectors in the plane, and choose 
a2 among the shortest lattice vectors in P not coUinear 
with ai and making an acute angle ai2 with it; so we have 
ai < 02 and, as argued above, 7r/3 < ai2 < tt/2. Because 
P is of highest density, the largest distance among lattice 
planes can be measured perpendicular to it. Accordingly, 
9b* = &3- Choose as among the shortest lattice vectors 
not contained in P and making an acute angle with at 
least one of ai or a2. One of the angles, ais of ai and as 
or 023 of a2 and as, can indeed be obtuse. However, if it 
is obtuse, we replace B by B* in the line of reasoning and 
continue with three acute angles: this will not influence 
the validity of Eq. (15). We have to examine two cases. 
First, suppose that ai = rs- Then 



from which the result follows by letting 77 go to zero. 



r-Bte* = 



27rai|ai x a2| 



27rai 



las • (ai X a2)| assina 



(29) 
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where a is the angle of a.^ to P. To maximize rsqE' , we 
choose 03 = ai = rs and then a to be minimum. None 
of ai3 and 0:23 can be smaller than ai2, otherwise the 
density of the plane spanned by as with either ai or a2 
would be higher than that of P. Therefore, the smallest 
a can be attained if as is in the bisector plane of ai2 
and ai3 = 023 = 0112 = tt/S. But then 02 cannot be 
larger than a\, otherwise, again, P was not a plane of 
maximum density. Thus, we conclude that ai = 02 = 
03 = rs and ai2 = 0:23 = 0131 = ""/S, specifying the fee 
lattice. The other case is tb = < ai. Then again 
rBQB* = 27r/sinQ:, but now a cannot be as small as 
before, otherwise the density of a lattice plane containing 



a3 would be larger than that of P. Thus, the maximum 
of rBQB' is indeed attained on the fcc-bcc pair. Its value, 
VGw, is easy to compute. 



C. Proof of the theorem, (i) 

Let X = u/^i(Sj + Yj), Xf C X finite, and let R be 
any finite configuration. If ip is continuous at each point 
of every B* then, making use of the definitions (l)-(3), 
(9) and the lemma. 



u{R\x \ Xf) = NRimpix) - mm + J m 



reR 



dk 



2(27r)'^ 



j=l OjtKeB^ reR 

Subtracting the corresponding expression in which Xf replaces R, we find 

U{R\X \ Xf) - txNn - U{Xf\X \ Xf) + nNx, = {Nx, - Nr) [h + ^(0)/2 - <^(0)p(X)] 

dk 



/ 



reR 



E' 

^eXf 



2(277)" 



OjiKeB' 



i,reR 



x-eXf 



r 



(30) 



(31) 



Suppose now that each qb* > Kq, where equality is al- 
lowed only if (,3 (k) = at |k| = i^o- Any such set of Bj s 
satisfies the continuity assumption. Indeed, if is continu- 
ous at K = and, because |K| > ifo for every nonzero K 
in yJ'j^iB* , If is continuous and takes on zero in all these 
points. Equation (31) is, therefore, valid, and the last 
term on its right-hand side vanishes for all R. The first 
term can be made zero cither by choosing Nr = Nxj or 
by setting ji = (p(O)p(X) — (p(0)/2. Because 1^ > 0, these 
choices prove that X is indeed a GSC and a /xGSC of ip. 

Next, we compute the energy density e{X) of X with 
each qb* > Kq. From Eq. (7) and the Lemma, 



Furthermore, if Ql denotes the cube of side L centered 
at the origin, then 



hm L-'' V 1 = p{X). 
xexnQL 



(33) 



J(x,X)=p(X)^(0). 



(32) With Eqs. (32) and (33), 



e(X) = 1 lim L-'^ ^ ^ p{^-^') = l ]im L-" ^ [/(x, X) - ^(0)] 

xexnQLX^x.'ex xexnQL 

= -Ip{X)p{0) + lp{X)'m = e{p{X)). (34) 

Here the first two equalities define the energy density of a general configuration, only the third one is specific to a 
GSC. We now prove that e(/o) is the minimum of e{X) among unions of periodic configurations of density p. As 
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earlier, we write an arbitrary (periodic or aperiodic) union of periodic configurations in the form X = Uj^i{Bj +yj). 
Computing I{x,X) with the help of Eq. (19), inserting it into the definition (34) of e{X) and separating the K = 
term we find 



e(X) = e(p(X)) + i ^p(i?,) ^(K)e'^ (''^-^-) hm L-'' ^ e^^'^. 



The limit can be evaluated: it yields p{Bj ) if K e B* n B* and zero if K e B* \ B* . Hence, we obtain 

J 

2 

2 



e{X) = e{p{X)) + 1^2 P(Bi)p{Bj) ^ <^(K)e'^-('^^-''' 

o#Kes?ns* 



OT^KeUBf 



^XB;(K)p(S,)e^^-^ 



(35) 



(36) 



where xb* (K) = 1 if K is in B* and is zero otherwise. 
Thus, e{X) > e{p{X)), as claimed. 

D. Proof of the theorem, (ii): absence of 
metsistability 

The following proposition makes no use of Fourier 

transforms, but exploits directly the summability as- 
sumptions (4) and (10). 

PROPOSITION, (i) Let be a strongly tempered pair 
potential. Let X and Y be locally uniformly finite config- 
urations such that the limits (33) and (34), defining the 
density and the energy density, exist. If p{X) = p{Y) but 
e{X) > e{Y), then X is not a GSC of (p, that is, X does 
not satisfy the local stability condition (5) restricted to 
number-preserving perturbations, (ii) If ip is not strongly 
tempered but summable on Bravais lattices in the sense 
of Eq. (10), and X and Y are unions of periodic config- 



urations with p{X) = p{Y) and e{X) > e(F), then X is 
not a GSC of ip. 

Proof. — The proof was already outlined in [1]. Here 
we give the missing details, (i) In the case of strongly 
tempered interactions consider the cubes Ql introduced 
above. In general, Al = NxhQl - NyhQl ^ but, 
because p{X) = p{Y), Al = o{L'^). Let Yl be a configu- 
ration in Qz, obtained from YhQl by adding or deleting 
\Al\ points so that Ny^^ = NxhQl- From the definitions 
(3) and (34) we deduce 

U{XnQL\X\QL) = e{X)L'^+h{XnQL,X\QL)+o{L'^)- 

(37) 

On the other hand, by strong temperedness 

U{Yl\X \ Ql) = U{Y n Ql\X \ Ql) + o{L'^), 
and therefore 



U{Yl\X \ Ql) = e{Y)L'' + I{Y n Ql, X \ Ql) - \l{Y f^QL,Y\ Ql) + o{L'^). 



(38) 



To conclude that 



U{Yl\X \ Ql) - U{X n Ql\X \ Ql) x [e{Y) - e{X)]L'^ < 



for L large enough, we have to show that the / terms in Eqs. (37) and (38) are of smaller order, o{L'^). Again, this 
holds because is strongly tempered. For example. 



L-''\I{Xr\QL,X\QL)\ = L-''\I{X nQ^_^,X\QL) + I{X nQL\QL-VL^X\QL)\ 



< p{x) 



sup Y k(x-r)| -hi ''/^sup ^ |v?(x-r) 

xeX,|x-r|>\/L xex 



+ 0(1) (39) 



r 



which indeed tends to zero as L goes to infinity, (ii) In the case when ip is not strongly tempered, let 
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X = u/^i(Bj + y,) and Y = i^jU{B'^ + y^-), with 
primitive vectors a.j^a and Sjo,, respectively. Con- 
sider the set U^=i{X^Q, '^aaj.a + Yjl I'^al < No], and 
let Xj be the union of J' non-overlapping adjacent 
translates of this set. Similarly, let R be the union 
of J non-overlapping adjacent translates of the set 
U/li{Ea'^"a;-„ + y^l \n^\ < 7V„}. In this way Nr = 
Nxf Defining Vr and Vx^ by the equalities U{R) = 
e{Y)VR and U{Xf) = e{X)Vxf, Vr/Vx,- ^ 1 as all 
Na 00, because p{X) = p{Y). On the other hand, 
from the convergence (10) it follows that the interac- 
tion terms are of smaller order, V^^I{Xf,X \ Xf) and 
V^^I{R,X \ Xf) tend to zero. This ends the proof of 
the proposition. 

The proposition implies the assertion (ii) of the The- 
orem, because any ground state configuration Y of den- 
sity p{Y) = p{X), described in (i), has an energy density 
e{Y) = eip{X)) < e{X). 

IV. SUMMARY 

In this paper we have described the ground state con- 
figurations of pair interactions (p and + tp, where ip is 
the inverse Fourier transform of a nonnegative function 
vanishing outside the sphere of radius Kq, and tp^sa, non- 
negative finite-range interaction of range tq. The GSCs 
of <p alone could be obtained, as already in [1], above 
a threshold density pd, while those of p + ip in a den- 
sity interval whose lower boundary is pd and the upper 
boundary is the close-packing density of hard balls of di- 
ameter tq. Below this density, tp alone would allow as a 



GSC any configuration with a nearest-neighbor distance 
not smaller than rg; combined with p, its role is to de- 
crease the degeneracy of the GSCs of ip. This reduced 
degeneracy is still continuous inside the interval, but at 
the boundaries p + ip has a unique GSC, which is, in 
three dimensions, the bcc lattice at the lower and the fee 
lattice at the upper density limit. This transition from a 
bcc ground state at p = ps to an fee one at p = \/2/ rg is 
the most interesting finding of the present work. 

The method used in [1] and in this paper cannot 
be applied to obtain the ground states of p below pd- 
The threshold value is a true critical density separat- 
ing the high-density continuously degenerate region from 
the low-density region in which the GSC is presum- 
ably unique, apart from Euclidean transformations. If 
1^(0) > 0, the analytic form of the relation between the 
density and the chemical potential also changes at pd [1] . 
We expect that at least in a subclass of interactions the 
unique bcc ground state at ps survives at lower densi- 
ties. Similarly, the mathematical method used here is 
not suitable to obtain the ground states of p + tp above 
the upper density limit. Again, this limit seems to be a 
true critical value, with a continuously degenerate ground 
state below and, probably, a unique ground state above 
it. This unique ground state may depend on the details 
of -0, but we expect it to be the fee lattice for a subclass 
of positive interactions. To clarify these questions, and 
also the nature of the curious liquid-like ground state be- 
tween the fee and bcc phases will be the subject of future 
research. 

This work was supported by OTKA Grants T 043494 
and T 046129. 
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